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We attempt to access the regime of strong coupling between charge carriers and transverse dynamics 
of an isolated conducting "stripe", such as those found in cuprate superconductors. A stripe is 
modeled as a partially doped domain wall in an antiferromagnet (AF), introduced in the context of 
two different models: the t-J model with strong Ising anisotropy, and the Hubbard model in the 
Hartree-Fock approximation. The domain walls with a given linear charge density are supported 
artificially by boundary conditions. In both models we find a regime of parameters where doped 
holes lose their spin and become holons (charge Q — 1, spin S3 = 0), which can move along the 
stripe without frustrating AF environment. One aspect in which the holons on the AF domain wall 
differ from those in an ordinary one-dimensional electron gas is their transverse degree of freedom: 
a mobile holon always resides on a transverse kink (or antikink) of the domain wall. This gives rise 
to two holon flavors and to a strong coupling between doped charges and transverse fluctuations of 
a stripe. 
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I. INTRODUCTION 

The nature of charge carriers in high-Tc cuprate super- 
conductors remains a subject of debate. The stoichiomet- 
ric ("parent") compounds are antiferromagnetic (AF) in- 
sulators well described by the Heisenbcrg model. Experi- 
ments show that, at least in some of the cuprates, doping 
with holes creates an intrinsically inhomogeneous state 
with periodically modulated charge and staggered spin 
densities. Neutron and X-ray scattering experiments in- 
dicate that staggered magnetization has a period of modsj 
ulation that is twice as long as that of charge density.lil'BI 
This is consistent with the notion of charged stripes sep^ 
rating AF domains with alternating Neel magnetization.cl 

Stripes as domain walls in the ground state of a 
collinearly ordered antiferromagnet were predicted — 
prior to reliable experimental detection — onJjhe basis of 
Hartree-Fock studies of the Hubbard models near half- 
filling. Mean-field calculations yield a linear density of 
1^=1 doped hole per lattice site, which almost cer- 
tainly means insulating stripes, in apparent contradic- 
tion with experiments. Besides, stripes observed in the 
cuprates tend to have a linear hole density «|J/2, 
at least when they are sufficiently well separated.cl To 
date, no reliable microscopic calculation yields the ex- 
perimentally observed filling fraction let alone explains 
the transport and high-temperature superconductivity in 
the cupratea__ Numerical simulations have so far been 
inconclusive .EfLl 

In the absence of a reliable microscopic theory, at- 
tempts have been made to find a phenomenological de- 
scription of the stripes. In cme of the more popu- 
lar routes, a stripe is modeledtl as a one-dimensional 
electron gaS|-|(lDEG) interacting with the surrounding 
enviropnieiitEl and with the transverse motion of the 
stripell30 As we argue below, this approach rests on 
the assumption that the physics of an isolated stripe is 
basically the same as in the limit z/ — > 1. In this limit, a 



stripe is almost completely filled with holes and can be 
described as an electron gas at low density 1 — i'. It is 
far from obvious, though certainly not implausible, that 
stripes with v w 1 and z^ w 1/2 should exhibit qualita- 
tively similar behavior. 

In this work we develop a qualitatively different (but 
not less plausible) phenomenology of a partially doped 
stripe. It is based on two different model calculations 
performed in the limit of low hole density on a stripe, 
I' ^ 1. Nominally, this is as far from the observed density 
1/ w 1/2 as the electron-gas limit v 1. The quantum 
numbers of charge carriers (holons) in our model calcu- 
lations are completely different from those of electrons. 
It likely means that the two limits are not adiabatically 
connected. It is clear then that the phase « 1/2 can 
resemble only one of the low-density limits: either — > 1, 
or z/ ^ — or possibly none of the above! 

Building on our model calculations we conjecture that 
charge carriers of the 'V — > 0" phase are holons (charge 
Q — 1, spin S3 ^ Q). In both models the loss of spin is 
compensated by the emergence of another spin-like de- 
gree pf| freedom, termed the transversal flavor by Zaanen 
et ali3 This happens because a holon always resides on 
a transverse kink or antikink of a domain wall. Thus 
holons are strongly coupled to transverse fluctuations of 
a stripe. Yet, the motion of such objects along the stripe 
is free, it does not produce any additional spin frustra- 
tion. Such a holon gas is clearly very different from the 
electron gas of the "v — > 1" phase. 

Of course, our approach should not be interpreted as 
a suggestion that a stripe with small linear hole density 
I' <C 1 can be stable in any model relevant to high- Tc ma- 
terials. On the contrary, a domain wall in an undoped 
antiferromagnet is a highly excited texture. A finite lin- 
ear density of holes is needed to stabilize a domain wall. 
In our analysis we always assume that the domain wall is 
supported externally {e.g., by the boundary conditions), 
while its untwisting is suppressed by a sufficiently strong 
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anisotropy (strictly linear polarization in our Hartree- 
Fock analysis). In a model where partially filled stripes 
appear in the ground state, neither assumption would be 
necessary. 

Indeed, we have reduced the symmetry of the problem 
in order to stabilize topologically the Ising-type domain 
walls.ll3 In any model with the 0(3) symmetry of the 
Neel order parameter, there can be no topological argu- 
ments for their stability. For example, domain walls can 
be continuously untwisted in a broad class of Ginzburg- 
Landau models with a continuous 0{N) symmetry; such 
domain walls are locally unstable. However, local insta- 
bilities are not an issue for globally stable configurations. 
In practical terms, untwisting does not occur if domain 
walls appear in the ground state of a system. In the 
context of high- Tc stripe phases, such nontopologically- 
stable domain walls were discussed in Ref. Their 
global stability requires frustration of the AF order on 
some microscopic or intermediate length scale, e.g., as a 
result of doping. 

Therefore, our model calculations should be viewed as 
an attempt to identify plausible ground states of an iso- 
lated stripe. In contrast to phenomenological approaches, 
we do not postulate effective one-dimensional models 
of a stripe. Instead, we derive them by starting with 
a two-dimensional model describing an antiferromagnet 
with a domain wall. Our 2D models may be unrealistic 
for the cuprates, but the resulting ID effective theories 
have the set of elementary excitations consistent with the 
paradigm of a stripe as a doped fluctuating domain wall 
in an antiferromagnet. In essence, we rely on universal- 
ity: if the number of qualitatively different ground states 
of a stripe (classified by_quantum numbers and spectrum 
of low-lying excitationsEj) is limited, all of them may be 
derived from simple 2D models. 

From this perspective, our work adds the ID holon gas 
to the list of potential stripe models. Note, however, that, 
in the presence of interactions, the description of a stripe 
in terms of holons, and a more conventional one in terms 
of electrons, are not necessarily incompatible. Both mod- 
els can be viewed as Luttinger liquids with different col- 
lective modes: charge and spin in the case q£, electrons, 
charge and transverse fluctuations for holons .Ej Thus, ID 
electrons with a spin gap and ID holons with a transverse 
gap may well represent one and the same phase. We in- 
tend to discuss the role of interactions among the holons 
in a future publication. 

The paper is organized as follows. In Sec. ^ we ana- 
lyze partially doped domain walls in a t-J model with 
large Ising anisotropy. The ground state and the spec- 
trum of elementary excitations (spinous and holons) are 



found explicitly in the limit of small doping ly. In Sec. Ill 
we present numerical evidence for holons in the Hubbard 
model, within the Hartree-Fock approach. The Hartree- 
Fock equations for the Hubbard model are further ana- 
lyzed in Sec. IV, where we introduce an appropriate long- 
wavelength approximation and study the spectrum of 
midgap states induced by a domain wall. We give heuris- 



tic arguments for the existence of fermion zero modes 
around a transverse kink on a domain wall. Technical 
details are collected in the Appendixes. 



II. HOLONS ON A DOMAIN WALL: T-J MODEL 
WITH ISING ANISOTROPY. 

A simple yet very instructive example of holon gas on 
a domain wall is offered by the t-~J model with king 
anisotropy, previously considered by Kivelson et al.pH 
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with the usual constraint of no double occupancies; the 
sum is taken over pairs of nearest-neighbor sites. The 
t-J model proper is restored if we set J± — Jz — —4 V. 

The analysis of the model (|^) is greatly simplified in 
the strongly-anisotropic limit Jz — > (X) (while the energy 
of an AF bond Ei, = V — Jz/4: may remain finite) . In this 
case, an isolated hole in the AF bulk is localized and has 
the energy e = —AEb, the cost of 4 missing bonds. Two 
holes on adjacent sites share one missing bond, which is 
to say that their interaction energy is Ei,. When severing 
a bond is costly {Ei, < and has a large absolute value), 
hole-rich islands are formed in the otherwise unaltered 
antiferromagnet, with the energy —2E[, per doped hole 
(phase separation in the bulk). In the opposite limit. 
Eh 3> doped holes strongly repel one another and 
stay apart. The energy per hole in this case ialZl 



ehoic = -4 Sfc = Jz - W. 



(2) 



Phase separation or not, doped holes are immobilized in 
the bulk of an antiferromagnet. As long as Jz greatly 
exceeds both t and J^, the cost of frustrated (ferromag- 
netic) bonds produced by a moving hole outweighs the 
reduction in kinetic energy. 

This situation changes dramatically in the presence of 
a domain wall: doped holes become mobile. 



A. Spinons on a domain wall 

Elementary excitations of an undoped domain wall are 
kinks formed in pairs by flipping two nearest-neighbor 
spms [Fig. |l| (a), (b)]. The kinks are mobile, carry zero 
charge and spin Sz = ±1/2. We thus tepBi them spinons. 
To make their relation to ID spinonsES more explicit, 
we have integrated spin Sz across the domain wall with a 
smooth envelope to obtain an effective ID spin chain rep- 
resenting the domain wall [Fig. |^ (a), (b), open symbols]. 
By using Bloch states,E3 one finds the spinon energy spec- 
trum, 
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i(fcx) = + cos2fc, + 0[Jl/J,) > 0. (3) 



Clearly, for ^ J±, such excitations are strongly 
gapped; an undoped domain wall is very stiff. 



(a) >/ A>/ A;^ A>/ A (c) AO A>/ A>/ A 



For small momentum fcj^, this is smaller then the energy 
of a single hole (||). Therefore, a dilute gas of holons has 
a lower energy than a collection of holes similar to that in 
Fig. |l| (c) . As shown in Appendix a dilute holon gas is 
stable against phase separation for Ef, > when the in- 
teraction between the holons is strictly repulsive. When 
Eb < 0, holes on adjacent sites attract. This attrac- 
tion wins over an increase in kinetic energy for Et ^ — 
causing phase separation: holes can lower their energy by 
forming a densely populated {v = 1) stripe leaving the 
rest of the domain wall undoped (see Appendix ^). 
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(b)>/ AA>/ A (d)^ AA>/ AO>/ A 

FIG. 1. Elementary excitations of a bond-centered domain 
wall in the anisotropic t-J model, p is the transverse flavor, 
(a-b) Two spin exchanges applied to the vacuum state create 
two separated spinons. (c-d) As a result of hopping, a doped 
hole leaves behind a spinon and becomes a holon. An anti- 
symmetric spin configuration around the holon means that it 
has = 0. 



B. Holons on a domain wall 

A single hole doped into a straight domain wall 
[Fig. 1^ (c), (d)] cannot hop along the stripe in the limit 
Jz oo: it must leave a spinon at the original location of 
the hole, which costs 0{Jz) in magnetic energy.c3 There- 
fore, a single hole can only oscillate across the stripe; the 
corresponding ground-state energy is 



Eho\c 



Jz 

^4Eb^f 



\t\ + o{tyjz). 



(4) 



It is important to realize that, once a hole starts mov- 
ing along the domain wall, it does not create a string of 
ferromagnetic bonds, which localize a hole in the bulk. A 
moving charge is now associated with a kink in the trans- 
verse position of the stripe [Fig. |l| (d)]. This composite 
object has spin Sz — and charge Q — +1. Following the 
spin-chain convention it can be termed a holon. More- 
over, holons can be created in pairs without spinons: two 
spinons created by two holes can annihilate each other. 
We can say that a localized doped hole decays virtually 
into a spinon and a holon. Another hole nearby can ab- 
sorb the costly spinon and become a holon. The interme- 
diate spinon is not needed if the two holes are on adjacent 
sites on the same side of the wall. 

A holon with momentum along the stripe has the 
energy 

£^hoio„(fcx) - -4^b - ^ - 2|t| cosfc, + 0{tyjz). (5) 



C. Preferred linear charge density 

So far we have considered an antiferromagnet with a 
single domain wall maintained by the appropriate bound- 
ary conditions. Within this model, one can study arbi- 
trary linear concentrations of holes v on the wall. Par- 
ticularly simple situations are the limit of dilute holes 
1/ ^ (gas of holons with a transverse flavor) and the 
opposite limit v 1 (ID electron gas). 

Because Jz is large, partially doped {u < 1) stripes 
do not occur naturally in this model. The preferred 
linear density of charge u aap be found by using the 
usual Maxwell construction.La To do so, one minimizes 
the energy per doped hole — including the cost of creat- 
ing domain walls. Since partially doped domain walls 
contain costly ferromagnetic bonds, they will not occur 
if ^ oo. A lower bound for the energy per doped hole 



e>{v-^-l)^-AEb + 0{t). 

In the limit we consider, this expression is a strictly de- 
creasing function of the optimal configuration of the 
stripe corresponds to v — \. 



D. Implications 

The anisotropic t-J model (^ dominated by the Ising 
term provides a good illustration to the strategy outlined 
in the Introduction. True, this simple model predicts in- 
sulating stripes with v = 1, contrary to experimental 
observations. Nevertheless, it has enabled us to find two 
possible phases of conducting stripes that may arise in 
more realistic models: the lILelectron gas (the 'V — + Vl, 
phase) suggested prcviouslyl3li3 and the ID holon gasll3 
(the 'V 0" phase). To do so, we have created a do- 
main wall by fixing boundary conditions and doped it to 
any given hole density v. At this stage, the simplicity 
of the model turns into a virtue: quantum numbers of 
elementary excitations can be readily determinedEj 

When macroscopic phase separation is absent, i.e., for 
Eb > [Jz < but t, J± <^ Jz), a gas of holons is 
formed on a weakly doped domain wall. Holons are mo- 
bile kinks of the domain wall with charge Q = +1 and 
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no spin, Sz = 0- Compared to their one-dimensional 
counterparts, domain- wall holons have an additional fla- 
vor, p — ±1/2 (isospin), which denotes the direction of 
the associated transverse kink [Fig. |] (d)]. An effective 
model describing the holon gas has been previously out- 
lined in Rcf. 

Even though it does not yield the observed stripe filling 
v ~ 1/2, the strongly anisotropic limit of the model (|l|) 
has certain appeal: it is simple enough to permit con- 
trolled calculations. In particular, spin waves are gapped 
and all associated dissipation effects are suppressed. In 
addition, holes are not allowed to leave the domain wall 
and therefore cannot go around each other. A domain 
wall in this limit is a strictly one-dimensional object, yet 
its transverse motion is fully accounted for. 

III. HOLONS ON A DOMAIN WALL: HUBBARD 
MODEL 

Could holons be generic to domain walls in an an- 
tiferromagnet or are they just a curiousity of the t-J 
model in the Ising limit? To answer this question, we 
have attempted to find similar excitations in the Hub- 
bard model. Clearly, the problem is much more difficult 
because there is no controlled approximation in this case, 
certainly not in two spatial dimensions. 

When the coupling strength U is weak compared to 
the free-electron bandwidth of At, one can hope to find 
some guidance in mean-field Hartree-Fock (HF) solu- 
tions. This approach has been successful — to a degree-pr 
in predicting the existence of stripes in the cupratesu. 
Numerical HF calculations show that, away from half- 
fiUing, doped charges form stripes along directions (0,1) 
or (1,1) with v = 1 doped hole per lattice site along a 
stripe. Such stripes are, indeed, AF domain walls. Be- 
cause they are filled with holes to capacity, charged ex- 
citations ^e gapped and thus an individual stripe is an 
insulator .E3 (The same problem arises in the t-J problem 
in the large- limit considered above.) 

Holons described in Section ^ are solitons with a well- 
localized charge distribution. To find analogous excita- 
tions in the Hubbard model (at the HF level) let us recall 
the properties of midgap states induced by domain walls. 
It is well known that solitons with anomalous quantum 
numbers arise in connection withfiermion zero modes in- 
duced by such topological defectsEil (see Appendix ^ . A 
uniform domain wall in 2 dimensions confines a midgap 
state only in one direction, across the wall. Solitons of fi- 
nite extension — both across and along the wall — can exist 
only if there is an inhomogeneity on the domain wall. In 
this section we shall show that a bond-centered domain 
wall with a wiggle may "bind" a holon, i.e., a soliton with 
quantized charge Q — 1 and spin S3 = 0. 

Although the solitons are static at the HF level, it is 
merely an artifact of the mean-field approximation, in 
which the average spin and charge densities are assumed 



to be time- independent. Mean- field configurations with 
solitons at different positions x along the wall should be 
viewed as degenerate minima of the action, i.e., as classi- 
cal solutions with a soliton at x. Quantization of the soli- 
ton restores broken translational symmetry: plane waves 
are superpositions of states with a soliton at all possible 
sites X. At-tlie semiclassical level, the energy of a soliton 
is given byEfl 

E = ^El+plvl^ »E,+ pIvI^/2Eo, (6) 

where Eq is the mean-field energy of the soliton and wid 
is the ID Fermi velocity calculated using the static HF 
wavef unctions. As discussed below, the holon energy 
spectrum is similar to that (H) of the t-Jz model, al- 
though both the inverse mass and velocity are substan- 
tially reduced (in the weak-coupling limit), reflecting the 
collective nature of the soliton. 

Just as in the case of the anisotropic t-J model, we 
are not dealing with the ground state of the model. A 
stripe needs a finite linear density of charge in order to 
be stable. Because collective excitations tend to be large 
at weak coupling, an appreciable linear density of holons 
likely requires the coupling to be strong, or else their 
overlap will completely destroy their individual proper- 
ties. With only a Hartree-Fock approach at our disposal, 
we cannot access the strong-coupling limit (although we 
have tried to mimick it in the anisotropic t-J problem). 
Instead, we maintain the domain wall by boundary condi- 
tions, and vary the linear charge density along the stripe 
by changing the total number of holes in the system. 

The weak-doping expansion in the Hubbard model, 
however, has an additional problem, which was not 
present in our analysis of the t—Jz model. Namely, 
the undoped domain walls are always unstable. Indeed, 
the undoped system can be accurately described by a 
Heisenberg-like model, and here the energy of a domain 
wall can be continuously lowered by perturbing in the di- 
rection orthogonal to the original magnetization vector. 
In the static HF configurations presented below [Figs. ^ 
^ , this untwisting instability was suppressed by imposing 
a constraint of linear polarization. Nevertheless, these 
textures are the (particular) solutions of the full set of 
HF equations. 

The domain wall (antiphase stripe) is favored by the 
holes. A fully-doped stripe at 1/ = 1 is both locally and 
globally stable already at the HF level. For a partially- 
doped stripe, this approximation (plus the constraint of 
linear polarization) gives static localized holes. We have 
found that the untwisting in the full set of HF equations 
(arbitrary polarization) starts to develop on the undoped 
portions of the stripe, compressing the remaining holes 
into segments of a fully doped stripe ending with semi- 
vortices similar to those discussed in Ref. We be- 
lieve that this is an artefact of the used weak-coupling 
approximation. Namely, we expect that at sufficiently 
strong coupling U the effective attraction between the 
holes (caused by the untwisting) will be compensated by 
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their increased mobility along the stripe.c3 In some sence 
this is similar to what happens in the anisotropic t-J 
model (|l]) in the region —t^Ei,< 0: even though holes 
can gain some potential energy by sitting next to each 
other, the associated loss of their kinetic energy prevents 
phase separation along the stripe. 

In the remainder of this section, we present numerical 
results obtained in HF calculations. Further analysis of 
the results in a long-wavelength approximation is given 
in the next section. 



A. Hubbard model: numerical results 

We have solved self-consistently HF equations of the 
Hubbard model 



remains equal to one everywhere, the soliton has zero 
charge (see Appendix ^). 



[a{v))-a-{p(v))]H^)=E4'(j), 



(the notation is explained in Sec. IV) for coUincar spin 
configurations. 



(si(r)} = (s2(r))=0, (S3(r))^0, 



(7) 



at small and intermediate interaction strengths U — 
2t . . . At. We always started with two AF domains hav- 
ing opposite values of staggered magnetization m. On 
the lattice row separating the two domains, staggered 
magnetization was initially disordered. Such configura- 
tions could thus later converge into site-centered, bond- 
centered, or meandering stripes. 



1. Undoped stripe 

At half-filling, the stripe always became bond-centered, 
as in the anisotropic t-J model. We have explicitly veri- 
fied that a site-centered stripe always has a higher energy 
in the absence of doped charges. 

In a few cases, an initially disordered stripe converged 
to a state with a higher energy, a bond-centered stripe 
with a defect where the domain wall shifts one lattice 
spacing sideways [Fig. §(a)]. Ripples in staggered magne- 
tization (— l)^"'"^(s3(r)) around the wiggle are an interfer- 
ence effect between staggered spin, varying as (—1)^+^, 
and a smooth compone^ of S3. By averaging {s^ir)) over 
four neighboring sites,EZl 

moo(r) = [(s3(r)) + {ssir + x)) + (s3(r + y)) 

H-(s3(r + x + y))]/4, (8) 

one can suppress the staggered component and uncover 
a spin soliton residing at the wiggle [Fig. ||(b)]. An even 
better view of the soliton is afforded when spurious long- 
range spin-density oscillations induced by the boundary 
are removed by plotting the symmetrized spin density 
(s3(r) -I- S3(— r)}/2 [Fig. 0(d)]. Because particle density 
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FIG. 2. A spinon on a bond-centered wall with a wiggle. 
CoUinear HF solution of the Hubbard model for U — 2t on a. 
square lattice 25 X 25. (a) Staggered magnetization m(r). (b) 
Smoothed spin density moo(r). (c) The deviation of total soli- 
ton spin S3{R), defined in Eq. (j^, from 1/2 for three different 
lattice sizes, (d) Symmetrized spin density (s3(r) + S3(— r)) /2. 

Numerically, the soliton has a total spin 5*3 = ±1/2: 



SsiR) = E E ^*3(r)) ^ ±1/2 as R 



(9) 



\x\<R\y\<R 



We have checked that the spin is well localized: S3{R) — 
53(00) vanishes exponentially with R [Fig. ||(c)]. This 
spin soliton is an exact analogue of the spinon found in 
the t-J problem [Fig. ||(b)]. Here we also find spinous 
of two flavors, those bound to transverse kinks and an- 
tikinks. 



2. Stripe with one doped hole 

With one hole added, an initially disordered stripe typ- 
ically converged to one of the two bond-centered config- 
urations. A straight bond-centered stripe would eventu- 
ally contain a polaron (Fig. ^), a nontopological soliton 
with the quantum numbers of a hole (Q = 1, 6*3 = ±1/2). 
The presence of a nonzero spin density is manifested 
in the typical ripples of the staggered magnetization 
(— l)^+^(s3(r)) — the result of the interference between 
the staggered and smooth spin components. 
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FIG. 3. A polaron on a straight bond-centered wall. 
CoUinear HF solution of the Hubbard model for U = 2t on 
a square lattice 25x25. (a) Staggered magnetization m(r). 
(b) Charge density 1 — n{r). (c) a;-staggered magnetization 
mio(r). (d) Smoothed charge density. 

The other typical configuration was a bond-centered 
wall with a transverse kink [Fig. ^, on which a charge 
Q = -t-1 is localized. The absence of interference fringes 
is a sign of zero spin density. Indeed, integrated spin 
Sz{R) is numerically zero (of order 10~^^) for all values 
of R. The absence of spin can also be verified by plotting 
77100 (r), the smoothed spin density (H). Since Q +1 
and 5*3 = 0, one immediately recognizes a holon in this 
soliton. As its counterpart in the anisotropic t-J problem 
[Fig. |^(d)], it binds to a transverse kink or antikink of the 
domain wall. 

Why does a wiggle on a domain wall change the nature 
of a doped charge in such a dramatic way? Essentially, 
a 6o77(i-centered domain wall can be thought of as a ID 
AF chain (roughly, two parallel spins across the domain 
wall create an excess spin 1/2). A straight domain wall 
corresponds to a spin chain with perfect AF order. If 
there is a transverse kink on the domain wall, the stag- 
gered magnetization of the effective chain changes sign 
at the transverse kink, i.e., staggered magnetization itself 
has a kink. In analogy 'wpth polyacetylene (as discussed, 
e.g., by Berciu and JohnEj), a doped charge becomes ei- 
ther a polaron (no AF kink), or a holon (an AF kink 
is present). To. illustrate this, we plot the x-staggercd 
magnetizationtJ 777io(r) 

777io(r)-(-l)-[(s3(r))-(s3(r + x)) 
+ (s3(r + y)) - (s3(r-|-x-Hy))]/4 

for a straight wall [Fig. |3|(c)] and for a wall with a wiggle 
[Fig. ^(c)]. In both cases, the x-staggered magnetization 
is confined to a narrow strip, which can be identified 
with the effective chain of the t-J problem. Clearly, 7tiio 
alters the sign at the center of a holon [Fig. ^(c)] but 



is only slightly depressed around a polaron [Fig. |^(c)]. 
Later, we will substantiate these qualitative arguments 
with an analysis of midgap states (particularly, fermion 
zero modes). 




FIG. 4. A holon on a bond-centered wall with a wiggle. 
CoUinear HF solution of the Hubbard model for U = 2t on 
a square lattice 25x25. (a) Staggered magnetization m(r). 
(b) Charge density 1 — n(r). (c) a;-staggered magnetization 
77iio(r). (d) The deviation of total charge Q{R) around the 
soliton from 1 for three different lattice sizes. 

Although the question of stability of large solitons on 
a weakly doped stripe is purely academic (see the dis- 
cussion in the Introduction), we have compared ener- 
gies of isolated polarons and holons. At weak coupling, 
U ^ 2.5i, polarons have a slightly lower energy. Two or 
more polarons preferred to bind into spinless bipolarons. 
For U > 3t, holons had a lower energy. Moreover, two 
well-separated holons (Fig. 5) had a smaller energy than 



a bipolaron for U > 3t. 



le trend is clearly to favor 



holons as the coupling gets stronger. 



p=-1/2 P=+"'/2 




FIG. 5. Staggered magnetization m{r) and hole density 
1 — 7i(r) in a HF calculation at U = 3t, 24 x 24 sites. A 
bond-centered wall with 2 wiggles and 2 doped holes. 
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IV. HUBBARD MODEL: A HARTREE-FOCK 
ANALYSIS 



-crit^V(r + Ar) - 0-3 ;7m(r) ■0(r) =E'iJj{r), (13) 

Ar 



To interpret the numerical results discussed in the pre- 
vious section, we have conducted a thorough analysis of 
electron states induced in the middle of the Hubbard gap 
by an AF domain wall. Midgap states of a straight do- 
main wall jhave been previously explained in great detail 
by Schulz.cJ Studying a domain wall with a wiggle on a 
lattice is a rather challenging task, therefore we first de- 
rive a long-wavelength approximation that could capture 
the essential physics (e.g., the difference between bond- 
centered and site-centered walls). We find that midgap 
states of an undoped domain wall resemble a half-filled 
ID chain with the Fermi momentum kp — Tr/2. The rel- 
evant low-energy states on a domain wall are composed 
of Bloch states with momenta near (±7r/2, ±7r/2). This 
leads to a theory of 8-component fermions (2 spin com- 
ponents X 4 Fermi points). Finally, we relate holons 
observed numerically to fermion zero modes induced by 
a wiggle, and trace the origin of holon transverse fiavor 
to the doubling of fermion components (8 instead of the 
usual 4 in the ID case). 



acquires a charge conjugation symmetry 

^(r) ^a2V*(r). (14) 

In this particle-hole symmetric form, the mean- 
field equatiopa r-resemble those in the theory of 
polyacetyleneB3'E3 The discrete symmetry (p^ ) has im- 
portant implications, among them the possibility of spin- 
charge separation. 

When staggered magnetization is uniform, the HF 
Hamiltonian (n3[) can be readily diagonalizcd using the 
momentum basis: 

El, = ±^el + U^m\ 

where Ek = — 2t(cos fc^; -l-cos fcj,). The one-electron energy 
spectrum has a gap 2Ao = 2J7|r7i|, which makes the sys- 
tem an insulator at half-filling. Note that, thanks to the 
spin axis rotation (|ll|), there is no doubling of the unit 
cell and the Brillouin zone is therefore not folded. 



A. Mean-field equations 



The HF equations for the Hubbard model are 
- i ^ V(r + Ar) - |[(a(r)) • a (p(r))] ^(r) = Ei^{v), 

(10) 



Ar 



where V'(r) is a 2-component spinor wavefunction, S is 
the triplet of Pauli matrices and the sum is over vectors 
Ar pointing to the four adjacent sites. The expecta- 
tion values of spin {<j{t)) and density (p(r)) should be 
calculated self-consistently. As discussed in the-. Intro- 
duction, we are specifically looking for collineaiBj'E3 so- 
lutions, therefore we set (cri) — {02) — 0. It is customary 
to rotate the spin axes on one of the sublattices through 
TT, which we do as follows: 



V'(r) 



x+y 



V'(r). 



(11) 



In the new basis, the HF equation reads 
- ait ^ ■0(1" + Ar) - U[asm{v) - n(r)/2] V(r) = E ipir) 

(12) 



Ar 



where mean- field parameters TO(r) and n{r) are the stag- 
gered spin and the charge density. 

To simplify further analysis, we neglect density 
fiuctuationf£il and set n{r) = n in Eq. (p^, also shift- 
ing E E + Un/2. The resulting HF equation. 



B. Straigiit domain wall along x 

A domain wall implies a change of sign for staggered 
magnetization m(r) along a line on the lattice. Because 
the "local gap" U\m{r)\ is reduced on this boundary, one 
expects to find electron states inside the Hubbard gap 
Aq. In the case of an isolated straight uniform domain 
wall, such states are localized across the wall and ex- 
tended along it. In view of the charge conjugation sym- 
metry ( [1^ ) , Eq. ([1^ ) has an even number of midgapj-sjtates 
for every momentum — tt < kx < tt, normally two.c3 



CM 



bond-centered 
site-centered 



**** 



^^^^ 









\(-H,+l) 




%^ 






(-1,-1).^ 


4+1-1) 



ji/2 



FIG. 6. Left: midgap one-particle spectrum E{kx) of an 
AF with a straight domain wall. Self-consistent solution of HF 
equations. U = 2.5t, 48x27 sites. Dashed line: low-energy ap- 
proximation, Eq. (po|). States outside the gap are not shown 
for the site-centered stripe. Right: 4 Fermi patches in the Bril- 
louin zone, Eq. (p^, shown with eigenvalues {t^,t^). Dashed 
line: Fermi surface of the noninteracting system. 
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Depending on the symmetry of the domain wall — it 
can be site or bond-centered — the two midgap bands in- 
tersect at kx — ±7r/2 or are separated by a smaller 
gap 2Ai, respectively (Fig. |^). Qualitatively, this can 
be understood by picturing electrons in the domain-wall 
bands as a ID electron gas in an external magnetic field 
Um{x,y) = -B(y) cos(7ra;), which is either antisymmetric 
(site-centered), or symmetric (bond-centered) in ?/. In 
the former case, B{y) averaged across the wall vanishes 
and we deal with effectively free ID electrons. In the 
latter case, the electrons feel a nonzero magnetic field 
staggered along x, which induces the smaller gap 2Ai 
near kx = ±7r/2. A more comprehensive discussion can 
be found in Appendix 

An undoped antiferromagnet is particle-hole symmet- 
ric, therefore only a half of the midgap states are filled. 
Since the total number of midgap states on a wall of 
length L is 2L, the maximum linear density of holes or 
electrons that can be doped is = 1. Once all the midgap 
states are filled, the system is again an insulator. 

Stability of the "magic" filling v = \ can be explained 
by the following qualitative argument. There is only one 
energy scale Aq in the limit of weak coupling. The pre- 
ferred filling is determined by minimizing the total energy 
cost — including the energy of domain walls — per doped 
particle. When v < \, holes are doped into a midgap 
band, where one-particle energies are much less than the 
Hubbard gap Aq — at least in the weak-coupling limit. 
Therefore, the main part of the energy cost comes from 
creating a domain wall, which should be of order Aq per 
unit length: 



E 



aAo 



V <1, 



a is of order 1 (in ID, a — 2/7r). Doping beyond v — 1 
puts holes into the lower Hubbard band, separated by 
the gap Aq, at a cost 



aAo + {v- l)Ao 



C. Continuum formulation 

In the weak-coupling limit, U <^ At, the character- 
istic distances are large, and a continuum approxima- 
tion of some sort should provide a sufficiently accurate 
description of the system. The approximation must be 
intelligent enough to tell apart, say, bond-centered and 
site-centered domain walls, which, as we have seen, have 
quite different one-particle spectra. The difference comes 
from a change in the symmetry of the domain wall, and it 
should be possible to describe within a continuum theory. 

To construct an effective continuum approximation for 
describing a weakly-deformed domain wall, we first have 
to identify the relevant modes. In a weakly-coupled Hub- 
bard model near half-filling, all low-lying excitations are 
concentrated near the Fermi- lines \kx ^ ky\ = tt . A do- 
main wall induces a ID midgap electron band, which is 
half-filled {kp = 7i'/2) if the domain wall is not doped. 
The relevant modes for describing a weakly-deformed 
domain wall should be located near the intersection of 
the two Fermi-lines, which gives four "Fermi points" 
k = (±7r/2, ±7r/2). Together with the spin index, this 
implies that the continuum description should be formu- 
lated in terms of 8-component Fermion wavefunctions. 
An alternative, more quantitative way to reach the same 
conclusion is presented in Appendix O, where a straight 
domain wall on the lattice is analyzed. 

We write an electron wavefunction as a sum of four 
terms with smoothly varying amplitudes: 



V's(l-) 



E E 

a=±l /3=±1 



V'a/3s(l")e 



i'K(ax-\-l3y) /2 



(15) 



1/ > 1. 



[As there is no folding of the Brillouin zone in our formal- 
ism, points (7r/2, 7r/2) and (— 7r/2, — 7r/2) are not equiv- 
alent.] In Eq. (p^), we have added two more indices, 
a — sgn kx and f3 — sgn ky to the staggered spin index s 
[Eq. (pi])] . Only those Fourier components of magnetiza- 
tion which connect the four Fermi patches are preserved: 



At I' = 1, the energy per doped hole has a cusp, where 
the derivative e'(j^) jumps from —aAo to (1 — a)Ao. The 
cusp is actually a minimum: since stripes with v = 1 are 
known to exist, they must have a lower energy per doped 
hole than the uniform AF state, e(l) < Aq, i.e., a < 1. 

The minimum of e(i/) may shift to a lower filling if 
there are two energy scales for carriers on a stripe or if 
a fairly large gap (comparable to Aq) opens up in the 
ID band at some value oi v < 1. An example of the 
former scenario, described by Nayak and Wilczek, gives 
a smooth minimum at an incommensurate filling deter- 
mined by the ratio of the two energy scales; this yields 
conducting stripes. In the latter large gap can 

possibly be the result of a commensurate filling z/q < 1, 
producing a cusp in e{i>); such stripes will likely be insu- 
lating. Neither argument provides a convincing explana- 
tion for experimentally observed conducting stripes with 
the "magic" filling fraction ly = 1/2. 



(S3(r)) 



m„;3(r)e*"(""+'3^). 

a=Q /3=0 



Each index has its own set of Pauli matrices, {rf }, {rf }, 
and {cr,;} for a, /3, and s, respectively {i — 1,2,3). Any 
two operators from different sets commute with each 
other because they act on different indices. Some of these 
operators have a transparent physical meaning: 

k«^(r3^r|), i-ir=rf, (-1)'^ = rf , (16) 



CTl 



S2 



172 e 



iQ r 



S3 



0-36 



iQr 



(17) 



Using this notation we write the continuum approxima- 
tion of the HF Hamiltonian (|l^) as 

Hhf = -2itaai{T^dx + T^dy) - Ua^imiy), (18) 
m = mil + "loirf + migrf + moorf rf. 
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where a is the lattice constant. Only TOii(r), staggered 
magnetization proper, exists in the bulk of the antiferro- 
magnet inducing the Hubbard gap A = [/|mii(oo)|. The 
energy spectrum near the 4 points (±7r/2, ±7r/2) has the 
form 



The other three components of m (e.g., the average spin 
density toqo) can be induced around defects only. As 
we will show, states localized on defects are particularly 
sensitive to these components. 



D. Midgap spectrum of a straight wall 

To warm up, let us derive the midgap spectrum of 
a straight domain wall along the x direction. Transla- 
tional invariance requires that toqi — rrioo = 0. From 
lattice solutions (Appendix |^) we know that the wall 
fermions have a gapless (gapped) energy spectrum for a 
site-centered (bond-centered) domain wall. The absence 
of a gap can be demonstrated by finding a fermion mode 
with zero energy at px = —idx = (i.e., kx = ±7r/2). 
The Schrodinger equation (nsh for i? = reads 



■ dy 



U 
2ta 



iTlmia{y)]ip{y). (19) 



If we neglect mio at first, solutions of Eq. (19) are eigen 



states of (72, Tg and, e.g., Tg (then each zero mode comes 
from a single Fermi patch). Eq. (^ reduces to 8 uncou- 
pled scalar equations, givingj-a total of 8 linearly inde- 
pendent solutions. As usual,EJ only half of them [those 
with eigenvalues cr2rJ'TOii(-|-cxD) < 0] are localized on the 
wall, so that there are 4 zero modes. 

Remarkably, in addition to the usual twofold spin de- 
generacy, there is another spin-like degree of freedom, 
which will prove to be the transverse flavor. The origin 
of isospin (at weak coupling) is thus exposed: compared 
to a ID chain, there are twice as many "Fermi points" 
on a straight domain wall in 2D — see Fig. ^, right. 

The difference between midgap spectra of site-centered 
and bond-centered walls arises already in the first order 
in miQ. Eq. (19) has four zero modes if rtuo is an odd 
function of y, i.e., for a site-centered walLcJ Otherwise, 
the middle band is split by a small gap (Fig. ^, left). 



1. Site- centered domain wall 

In the presence of a nonvanishing mio, the spectrum 
E(px) of a straight site-centered stripe can be determined 
approximately by starting with Px = [Eq. (^9|)] and 
treating the term containing dx in the Hamiltonian ( p^ ) 
perturbatively. In the limit p^; — > 0, states outside the 
main gap can be neglected, which reduces the Hilbert 
space to the four zero modes [Eq. (f9|)]. By using the 



degenerate perturbation theory, we find a Dirac spectrum 
(dashed fines in Fig. ^, left): as Px — > 0, 



E{px) ±vpx, V = 2ta{Tf). 



(20) 



This compares well with a similar result (C7) obtained 
on a lattice. 



2. Bond-centered domain wall 

Alternatively, one can explore the limit of a small x- 
staggered magnetization, mio ^ mii(oo). In that case, 
by starting with mio = 0, one finds 4 degenerate zero 
modes for any px- A nonzero mio induces a splitting of 
the zero modes. To lowest order in mio, the energy at 

= is 

£;(0) = ±Ai - ±C/(mi) =±U [ mw{y)^\y)^{y)dy. 



As claimed, the gap is proportional to the x-staggered 
magnetization felt by an electron on the domain wall. 



E. Zero modes at a wiggle 

One way to prove that a charged soliton at a wiggle 
indeed has zero spin is to show that the HF equations 
contain a doubly degenerate fermion zero mode. An ele- 
mentary discussion of the connection between zero modes 
and separation of spin and charge is given in Appendix 
Because the problem is essentially two-dimensional (the 
domain wall is curved), it is much harder than its ID 
analogs. In ID, symmetry arguments arc generally suf- 
ficient to prove the existence of zero modes in ID — even 
on a lattice! (See Appendix |c[) In contrast, we have not 
been able to find such a general proof for the 2D problem 
of a wall with a wiggle — neither on the lattice, nor in the 
long-wavelength approximation. 

Instead, we offer a somewhat hand- waving argument in 
favor of zero modes in this case. Lack of rigor is compen- 
sated by an insight into the origin of the transverse flavor: 
it turns out that holons residing on transverse kinks and 
antikinks come from different points of the Brillouin zone, 
i.e., they are made of completely different stuff. 



> y X y > 

y A y > y 



y 



^ / ^ / ^ / ^ )^ 



FIG. 7. Left: Bond-centered stripe with a wiggle as a su- 
perposition of a site-centered domain wall (black arrows) and 
a ID AF chain with a kink (open arrows). Right: Holons with 
alternating isospins form a cite-centered stripe. 
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As illustrated in Fig. magnetization on a bond- 
centered wall with a wiggle can be obtained by super- 
imposing m(r) of a straight site-centered wall and that 
of a spin chain with a kink in x-staggered magnetiza- 
tion. Away from the wiggle, moo(r) — ™oi(r) — 0. To 
simplify the discussion, we will neglect these components 
altogether (but this is nonessential and can be remedied). 
Decompose m(r) into an a;-independent part and the rest: 

m(r) = m'-^\y) + m'-^\r), 
The Hamiltonian (|8|) can now be split in two parts: 



(0) 



jjp - -~2itaaiTldy -Uazlm-^^iy) -f jn-l'^ {y)Tl], (21) 
H^^l = -2ita - Ua^ [m ^ (r) + m'^^^ (r)rf ] . (22) 

As shown above, the "transverse part" (|2^) has 4 zero 

modes for each p^. Within this Hilbert space, iJ^F 
scribes right and left-moving fermions with spin, which 
see a staggered magnetization 

(mi(.T)) = j dy v){y)[mia{Y) + mii{Y)Tf]u{y) , 

where u{y) is a zero mode (^ of Eq. (|l|). The 
midgap fermion band acquires a gap of its own, Ai — 
?7|(mi(oo))| < A, with two zero modes (one for each 
spin) inside this smaller gap. "Longitudinal" wavefunc- 
tions of the two zero modes satisfy the equation 



,(0), 



(72- 



dip{x) 



U 



dx 



2ta{Tf) 



(23) 



The existence of two holon flavors can now be deduced 
from Eqns. (19) and ([2^). The zero modes have a finite 
norm only if 

aarfm^^i^ (+oo) < 0, (T2T3"(mi(+oo))/(Tf ) < 0. 

It follows then that the product of eigenvalues 

Tg^rl = sgn[TO(")(y = +oo) (mi (a; = +oo)) (rf )] (24) 

can be identified with the holon isospin 2p. This can 
be seen by extrapolating Eq. (^) to larger values of 
U, which reduces the size of holons. We have (rf) = 
((— 1)^) = (—1)^", where ?/o is the row number of the 
chain in Fig. |^. According to Eq. (p^, if rf = -1-1, 
spins on the chain and to the right (left) of the wiggle 
are an extension of the upper (lower) AF domain, as for 
the p = +1/2 wiggle in Fig. 0. Thus, p = t§t^/2. This 
identification is consistent with numerical HE solutions 
(Fig. ^), where rf Tg — sgnk^ sgiiky can be inferred from 
the orientation of a holon — the perfect nesting of the 
Fermi surface makes holon wavefunctions cigar-shaped 
and oriented along a lattice diagonal. 



V. CONCLUSIONS 

We have attempted to infer a set of plausible quantum 
numbers of low-lying excitations on a partially doped do- 
main wall in a strongly-correlated antiferromagnet by an- 
alyzing artificially-created domain walls in simpler sys- 
tems. Specifically, we have studied quantum numbers 
of well-separated holes doped inta domain walls in the 
t-J model with Ising anisotropyEZi and in the Hubbard 
model (in a Hartree-Eock approximation) . In addition to 
a usual ID electron gas,Q we have identified a new poten- 
tial candidate: the ID gas of holons. In this phase, which 
we have found at sufficiently small linear hole density z/ 
in both models, charge carriers (holons) have spin S3 = 
and charge Q — +1. Each holon resides on a transverse 
kink of the domain wall, which leads to a strong interplay 
between charges and transverse fluctuations of a stripe. 

We find it very encouraging that the charge carri- 
ers with identical quantum numbers result from two 
vastly different calculations. In the strongly coupled t- 
Jz model, the holons are small and immediately evident 
[Fig. 0(d)]. In the weakly-coupled Hubbard model, they 
are large and represent a collective effect (fermion zero 
modes). This indicates that a universality of some sort 
is at play, and, therefore, that the same "i/ — > 0" phase 
could result from more authentic models. Whether or not 
this phase is relevant for the cuprate stripes, which have 
v w 1/2, remains an open question. In future, we intend 
to extend this work to intermediate linear hole densities. 
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APPENDIX A: MORE ON THE T-Jz MODEL 

In this appendix we give some estimates of the en- 
ergies of domain walls in the strongly anisotropic t-J 
model (^. In particular, we show that in this limit doped 
holes can lower their energy by forming (fully-packed) do- 
main walls. This possibility has not been considered in 
Ref. |l^. We also discuss stability of the dilute holon gas 
considered in Sec. ||. 

In the limit of infinite J^, only fully doped domain 
walls, similar to those considered by Osnian et al.l^ have 
a finite energy (all frustrated bonds are covered with 
holes). A domain wall must therefore maintain continu- 
ity, i.e., adjacent holes must be nearest or next-nearest 
neighbors. This implies that a domain wall horizontal on 
average, can change its height y{x) by at most one unit 
at a time. Such a domain wall can be fully described 
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using a one-dimensional language, namely by specifying 
differences in the heights of neighboring holes: 

y{x+ 1) ~y{x) = -1,0, or 1. 

The system is therefore equivalent to a spin-1 chain ,@ 
with the Hamiltonian 



H^s - 



-3 E, - E, + - {St S-^, + h.c.) 



(Al) 



where n — x + 1/2 and E\j — V — Jz/4: is the energy 
of an AF bond. The first term in Eq. (Al) represents 
the energy of a straight segment of a wall (three broken 
bonds per hole), the second term counts the number of 
additional broken bonds due to kinks, while the last term 
describes the transverse hops of holes. 

The properties of the effective spin Hamiltonian 
and its generalizations have been extensively studied.SS 
When Eb » t, the AF ordering = (-1)" (a zigzag 
wall) is favored. For Ei, large and negative the ground 
state corresponds to S^-^ = (a flat domain wall). At 
\Eb\ ^ the system enters an intermediate critical phase, 
in which the density of domain wall kinks varies contin- 
uously. 

Up to terms of higher order in t/\Eb\, the energies (per 
hole) of the two ordered phases are 



£fiat — —3Eb — /\Eb 



^AEb-t'/Eb. (A2) 



For Eb < (two holes attract), phase separation in the 
bulk affords a lower energy, —2Eb < £flat, thus hinder- 
ing natural formation of stripes. For a strongly repul- 
sive bond, Eb 3> stripes win over a lump of immo- 
bilized holes: transverse fluctuations of a zigzag stripe 
lower its kinetic energy by an amount of order t per hole 
[cf. EcL (||)]. This possibility has not been considered in 
Ref. |l|. 

Let us now consider a domain wall created artificially 
(e.g., by boundary conditions at infinity) in order to 
study partially doped domain walls. In the absence of 
holes, such a wall is &ond-centered and straight to min- 
imize the number of broken bonds [see Fig. The 
corresponding energy cost per unit length is Jz/2. 

When holes are added to an antiferromagnet with a 
domain wall, they will necessarily bind to it (each hole on 
the domain wall reduces the number of frustrated spins; 
the energy gain is ^ Jz/2 per hole.) As discussed in 
Sec. |l[ a single doped hole acquires mobility by riding 
a kink [See Fig. |l|d]. The corresponding energy is given 
by Eq. (||). Assuming that holons are well separated, the 
energy per added charge is 



£holon gas = -4.Eb -2t + 0{f/Jz, V t) . 



(A3) 



Here we have not included the energy cost of creating a 
domain wall. 



The assumption of large separation may be violated 
even at very small linear densities i/ <C 1 if there is an 
attractive interaction between holes. For example, for 



Eb ^ —t, the energy (A2) of the fully-packed flat phase 



can be smaller than that of the holon gas, ( AS ) . The holes 
on the stripe will separate into a dense phase (y = 1) 
leaving a portion of the domain wall comletely undoped 
(^ = 0). 

To find a strict upper bound E'™™, such that for Eb < 
i?™'" the phase separation definitely happens, we can use 
a variational estimate for the ground state energy of the 
Hamiltonian (Al). The simplest estimate corresponds to 
all — 0, which immediately gives e < —3Eb. This is 



smaller then the energy of a dilute holon gas (A3) for 



Eb < i?™™ = —2t. This estimate of the phase separation 



boundary ii^™™ can be easily improved (increased) by 
using more sophisticated variational wavefunctions. 

On the other hand, phase separation of this sort is not 
expected for Eb > 0, when holes have uniformly repulsive 
interactions. This statement can be made more formal 
by evaluating a strict lower bound on the energy of any 
dense hole phase described by the Hamiltonian (Al). Es- 



timating each term in the Hamiltonian independently, we 
have, for Eb > 0, 



■AEb — 2\t\) = — A'';! Eholon gas- 



The inequality is strict because the terms in the original 
Hamiltonian do not commute. It implies that a dilute 
holon gas is stable to phase separation to a completely 
doped region and a region of an undoped stripe, as ex- 
pected on physical grounds for a repulsive interaction. 



APPENDIX B: ZERO MODES AND 
SEPARATION OF SPIN AND CHARGE 

We retrace the relation betweea-iurmion zero modes 
and separation of spin and chargeH3l23 For completeness, 
we will use the lattice version of the HE equations (|lj) , 

-crit^7/'(r + Ar) - (73 C/m(r) V(r) = E%lj{y). (Bl) 



1. Symmetries 

So long as we deal with coUinear AF configura- 
tions, one component of spin — ((T3/2) exp (z Q-r) in our 
notation — is a conserved quantity. The transformation 



V'(r) 



0-3 e 



iQr 



^(r) 



(B2) 



is a symmetry of the mean-field Hamiltonian. Here Q = 
The unitary part of charge conjugation (1#, 
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■0(r) a2il^{r), 



(B3) 



alters the sign of E and as such is not a symmetry of the 
mean-field equations. Rather, it can be referred to as a 
symmetry of zero modes. 

For the sake of convenience, we also want the system 
to be reasonably symmetric with respect to some sort 
of a parity transformation r r. If, however, the 

system has a domain wall, staggered magnetization will 
be antisymmetric under parity. To make it a symmetry, 
we combine parity with a spin flip: 



(B4) 



This "combined parity" is a symmetry of HF equations, 
provided that m(— r) = — m(r). 



2. Fermion zero modes 



Under some circumstances, Eq. (Bl) has solutions with 
zero energy. Such solutions are normally localized on a 
topological defect, e.g., on a domain wall. Note that a 
straight domain wall confines a fermion zero mode in one 
direction only — across the wall. Therefore, a soliton of 
finite dimensions requires a domain wall with an inho- 
mogeneity of some sort, e.g., a wiggle. Here we will as- 
sume that the wavefunction of a zero mode decays quickly 
enough in all directions. 

Zero fermion modes always come in doublets. This is 
essentia lly a consequence of the charge conjugation sym- 
metry (B3) at i? = 0. More formally, by starting with 
the symmetries of zero modes ( ]b^ ) and (B4), we can 
construct a triplet of SU(2) generators 



S2 



11 
2 ' 



iQ r 



Si = -i[S2, 5*3] = 



cTie 



(B5) 



Components of S are conserved quantities for a zero 
mode. By inspection, S-S — 3/4, i.e., a zero mode is 
a doublet (5* = 1/2). Physically, 6*3 is the component of 
the total spin of the system parallel to staggered magneti- 
zation and is therefore a conserved quantum number. Si 
and S2 are components of the total staggered spin. They 
generate staggered rotations of spins and are conserved 
for zero mode fermions, but not for bulk states. 



3. Spinons and holons 

We consider a system with only one zero-mode doublet 
localized on a topological defect at the origin, r = 0. 
We assume that the system is symmetric, i.e., to(— r) = 
— m(r). 

Half-filled system. First let us show that a half-filled 
system has a uniform density. States with E < Q are 
filled, while those with E > Q are empty. One of the two 



zero modes is filled. Because the density contribution 
of the zero mode ■!/'g(r)?/;o(r) is invariant under rotations 
generated by S, we can consider the situation when the 
mode with 5*3 = 1/2 is occupied without loss of general- 
ity. Then the operator 02 toggles between occupied and 
unoccupied states. The expectation value for the density 
is 

occ occ 

n(r) =^^t(r)^(r) =^V^(r)a2V(r) 



unocc 



all 



Thus a half-filled system has a uniform charge density. It 
will be shown below that it contains a charge-0, spin-1/2 
soliton at its center. 

Half- filled system ±1 electron. When a single electron 
or hole is added, the density distribution n(r) is deter- 
mined by the profile of the zero mode V'o(r)V'o(r). Be- 
cause the wavefunction is localized around a defect at the 
center, we find a soliton with charge ±1. 

At the same time, the soliton has zero net spin 5*3. 
Moreover, 6*3 = in any symmetric finite area around 
the defect. This happens because contributions to the 
total spin from r and — r cancel each other: 

(s3(-r)) = (-ir+^m(-r)= -(-l)--+^m(r) = -(,S3(r)). 

Thus, at half-filling ±1 electron, the system has a soli- 
ton with charge Q = =Fl and spin 5*3 — 0. Exactly at 
half- filling, the soliton has Q — Q and S^ = ±1/2. These 
are, respectively, a holon and a spinon. 



APPENDIX C: STRAIGHT DOMAIN WALL: A 
LATTICE ANALYSIS 

Consider a straight domain wall along the x axis. At a 
given lattice momentum along the wall, the mean-field 
Hamiltonian (O) reads 



ait[i>{y + 1) + iP{y - 1) + 2cosKi:{y)] 
-a3Um{y)tp{y) = Eilj{y). 



(CI) 



For definiteness, the domain wall is located on the line 
2/ = 0, so that m{--y) = —m{y). The site indices are in- 
teger (0, ±1, . . .) for a site-centered wall and half-integer 
(±1/2, ±3/2, . . .) for a 6o7irf-centered wall. 

We will first show that the smaller gap (separating the 
midgap bands) is absent in the case of a srfe-centered 
domain wall, for which m(— j/) = —m{y). It suffices to 
show that there exists a zero mode at = ±7r/2. 



1. Site-centered wall: zero modes 

Setting E — and = ±7r/2 converts Eq. (Oh into 
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+ 1) + ij{y - 1) = i(T2A(y) 
where A(j/) = Um{y)/t. In the bulk, 

A(+oo) = -A(-cx)) = Ao > 0. 



(C2) 



We are looking for a finite-norm solution to Eq. (C2). 
Evidently, a2 can be immediately diagonalized: (T2 — il- 
The substitution 

where x is an arbitrary constant spinor, yields a differ- 
ence equation with real coefficients for a scalar (f>{y)'. 



cj){y + l)-cf>{y-l) = -A{y)cj){y). 



(C3) 



Eq. ( p3D has two linearly independent solutions with 
the asymptotic behaviors 

4'±{y) ^ exp (±Ky) as y ^ +oo. 

The spread of the wavefunction is given by the equa- 
tion 

2 sinh K = Aq. 

Clearly only (y) can be normalizable — provided that 
it is also well behaved as y ^ — oo. It will be seen shortly 
that this is indeed the case. 

Observe that one can write 0_ (y) as a superposition of 
another pair of linearly independent solutions, the eigen- 
statcs of parity, 

where i — 1,2 and rji ~ ±1. Ordinarily, 771 ^ rj2, so that 

<f>-{y) = ci0i(y) +C202(y) 

is not, in general, a parity eigenstate. Yet it must be 
if it is to have a finite norm. Indeed, from different 
asymptotic forms of the solutions 4>±{y) we infer that 
(j)-[—y) = rj<j)-{y), or else it will diverge at y = —00. 
Luckily, all solutions of Eq. (C3) have the same (even) 
parity for a siie-centered domain wall: A(0) = and 
therefore 



(1) 



^1) 



-A(0)(?!)(0) = 0. 



Since 4>-{y) vanishes exponentially at both y — +00 and 
— cx), it has a finite norm. 



We have thus proven that the HP Hamiltonian (CI) 
has 2 solutions of a finite norm with = for \kx \ = 7r/2, 
one for each eigenvalue of (T2, if the wall is s?ie-centered. 
As a rule, there are no zero modes if the domain wall is 
6oTi(i-centered. 

The asymptotic behavior of the zero modes as y — > 00 

is 



'0(r) ~ exp 



X - « 2^22/ 



where 2 sinh k = Aq and x is an arbitrary constant 
spinor. We stress that, in the limit of weak coupling 
K <C 1, the Fourier components of V'(r) come from the 4 
"Fermi patches" \kx\ ~ \ky \ w 7r/2. This means twice as 
many Fermi points as in an ordinary ID electron gas! It 
is for this reason that holons (and spinous) on a domain 
wall come in two flavors. 



2. Site-centered wall: ID electron band 

For kx / 7r/2, the degeneracy of the two midgap states 
is lifted by the additional term —2taiC0skx in the HF 
Hamiltonian (CI). Near k^ — ±7r/2, this term can be 
treated as a perturbation acting in the Hilbert space of 
the two zero modes, allowing us to determine their split- 
ting to first order. 

Parametrize the wavefunctions in the familiar way. 



'ip{x,y) 



>(2/)' 



,-i(7r/2)o-2y 



X, 



where 0(y) is the solution of Eq. (C3) with the norm 1. 
The resulting i}}{x,y) diagonalizes the mean- field Hamil- 
tonian ( |Cl| ) with coskx = 0. The two-fold degeneracy 
of the zero mode is due to the freedom in choosing the 
spinor x- the framework of the degenerate perturba- 
tion theory, we compute the matrix element of the per- 
turbation = — 2<CTi cos/cj; between the states \kx\, Xi) 
and \kxi,Xi)'- 

'\kxi,Xi\ii\\kx\,X\) = -2tcoskxi 

X 2nS{kxi ~ kx2) (xkiXi) Y.^-lY<f'{y). (C5) 

y 

This result is quite tangible: as far as the midgap states 
are concerned, "integrating out" the transversal degree of 
freedom y yields electrons with nearest-neighbor hopping 
along the domain wall and no staggered magnetization — 
Cf. Eq. (^: 

- aitiY)[il,{x + 1) + ij{x - 1)] = E^{x). (C6) 
The effective ID hopping amplitude is 

t,^ = t{{-i)y)^tY,{-^?^\y)- (C7) 

y 

Note that tiD ^ i in the limit of weak coupling. 
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